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I. INTRODUCTION: EVENTS OCCUR IN 
COMPLEX TIME 



Instantons generally refer to solutions of classical equa- 
tions in the Euclidean space-time, i.e. once a Wick ro- 
tation t i y —it has been performed on time t in the 
Minkowskian space-time. Since the mid-seventies, they 
have been extensively used in gauge field theories to de- 
scribe tunnelling between degenerate vacua [l| . In intro- 
ductory texts [2|, y, for instance] , they are first presented 
within the framework of quantum mechanics: when the 
classical Hamiltonian of a system with one degree of free- 
dom has the usual form 



H (p, q) = \v 2 



V(q) 



(1) 



the original (non Wick-rotated) Hamilton's equations — 
may manifest itself through a transition that is not ne- 
cessarily a classically forbidden jump in position [J]. For 
instance, the reflection above an energy barrier, as a for- 
bidden jump in momentum, is indeed a tunnelling pro- 
cess [Hill]. In the following, we will consider the case 
of a simple pendulum whose dynamics is governed by 
the potential V(q) — 7COsg; for energies larger than 
the strength 7 > of the potential, one can observe a 
quantum transition between states rotating in opposite 
direction while two distinct rotational classical solutions, 
obtained one from each other by the reflection symmetry, 
are always disconnected in real phase-space. 

A second example is provided by a typical situation of 
resonant tunnelling: when, for instance, V has a third, 
deeper, well which lies in between two symmetric wells 
(see Fig.EK), the oscillation frequency between the latter 
can be affected by several orders of magnitude, when two 
eigenenergies get nearly degenerate with a third one, cor- 
responding to a state localised in the central well. Then, 
we lose the customary exponential weakness of tunnelling 
and it is worth stressing, coming back for a second to 
quantum field theory, that having a nearly full tunnelling 
transmission through a double barrier may have drastic 
consequences in some cosmological models. In this exam- 
ple, one can immediately see (Fig. [3b) that —V also has 
an energy barrier and working with a complete Wick ro- 
tation only remains insufficient in that case. 

In order to describe the tunnelling transmission at an 
energy E below the top of an energy barrier, which may 
be crucial in some chemical reactions, the pioneer works 
by Freed Q , George and Miller [1, 0] , have shown that 
the computation of a Green's function G(qf, \ E) (or a 
scattering matrix element) requires taking into consider- 

ation classical trajectories with a complex time. These 

complex times come out when looking for the saddle- 
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(p and q denote the canonically conjugate variables), the 
Wick rotation induces an inversion of the potential and, 
then, some classical real solutions driven by the trans- 
formed Hamiltonian p 2 /2 — V(q) can be exploited to 
quantitatively describe a tunnelling transition. As far 
as we know, in this context, only the simplest situa- 
tions have been considered, namely the tunnelling de- 
cay from an isolated minimum of V to a continuum and 
the tunnelling oscillations between N degenerated mi- 
nima of V that are related by an iV-fold symmetry In 
those cases, what can be captured is tunnelling at the 
lowest energy only. However, not to speak of the highly 
non-trivial cases of tunnelling in non-autonomous and/or 
non-separable multidimensional systems, there are many 
situations that cannot be straightforwardly treated with 
a simple inversion of the Id, time-independent, potential. 

First, tunnelling — i.e. any quantum phenomenon 
that cannot be described by real classical solutions of 
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time propagator, 

G(q f , qi ;E) = -= G(q f , qi ;t)e iEt ^dt, (2) 
V ZttTl Jo 

which is, up to now, the common step shared by all 
the approaches involving complex time [IoHl4l for ins- 
tance]. Though well-suited for the study of scattering, 
indirect computations are required to extract from the 
poles of the energy Green's function © some spectral 
signatures of tunnelling in bounded systems. Generically, 
these signatures appear as small splittings between two 
quasi-degenerate energy levels and can be seen as a nar- 
row avoided crossing of the two levels when a classical 
parameter is varied [l5l - [l8j |. In the present article, we 
propose a unified treatment that provides a direct com- 
putation of these splittings (formulae [9] or [10]) ; as shown 
in section [Til it takes full advantage of the possibility of 
working not necessarily with purely imaginary time, but 
with a general parametrisation of complex time as first 
suggested in [Tjj. The semiclassical approach naturally 
follows ( sec. IIIip and some general asymptotic expansions 
can be written (equations (IC6[) and (|¥U|)) and simplified 
(equation (17310 : they constitute the main results of this 
paper. To understand where these formulae come from 
and how they work, we will start with the paradigmatic 
case of the double-well potential (sec. IIVI) and the sim- 
ple pendulum (secE}; we defer some general and tech- 
nical justifications in the appendices. Then we will treat 
the resonant case in detail in sec. I VII where an appro- 
priate incomplete Wick rotation t i-> cr l6 t provides the 
key to showing how interference effects a la Fabry-Perot 
between several complex trajectories reproduce the non- 
exponential behaviour of resonant tunnelling, already at 
work in open systems with a double barrier [20j, |21| . After 
having shown how to adapt our method to the computa- 
tions of escape rates from a stable island in phase-space 
(sec. IVII[) . we will conclude with more long-term con- 
siderations by explaining how our approach provides a 
natural and new starting point for studying tunnelling in 
multidimensional systems. 



II. TUNNELLING SPLITTINGS 

A particularly simple signature of tunnelling can be 
identified when the Hamiltonian has a two fold sym- 
metry and, therefore, we will consider quantum systems 
whose time- independent Hamiltonian H = H(p, q) com- 
mute with an operator S such that S 2 = 1 (the " allows 
to distinguish the quantum operators from the classical 
phase-space functions or maps). In most cases, S stands 
for the parity operator: 

H(-p,-q) = H(p,q) . (3) 

The spectrum of H can be classified according to S and, 
for simplicity, we will always consider a bounded system 



whose discrete energy spectrum and the associated or- 
thonormal eigenbasis are defined by 

H\<fi)=E±\<fi) ; S|^}=±|0±) (4) 

where n is a natural integer. When the Planck cons- 
tant h is small compared to the typical classical actions, 
standard semi-classical analysis [22l424l | shows that one 
can associate some classical regions in phase-space to 
each eigenstate \<p^). This can be done by construc- 
ting a phase-space representation of |<^„), typically the 
Wigner or the Husimi representation, and look where the 
corresponding phase-space function <jx^ (p, q) is mainly 
localised (all the more sharply than H is small). For 
a Hamiltonian of the form (fTJ where V has local mi- 
nima, some of the eigenstates remain localised in the 
neighbourhood of the stable equilibrium points. For ins- 
tance, for a double-well potential whose shape is shown 
in Fig. d^), the symmetric state |</>+) and the antisym- 
metric state \4>n) a t energies below the local maximum 
of the barrier have their Husimi representation localised 
around both stable equilibrium points (p, q) = (0,±a), 
more precisely along the lines H(p,q) — E (Id tori) at 
energy E ~ £+ ~ E~. Because of the non-exact de- 
generacy of E^ and E~, any linear combination of 
and constructed in order to be localised in one well 
only will not be stationary anymore and will oscillate 
back and forth between —a and a at a tunnelling fre- 
quency \E~ — E^\/(2nh). At low energy this process 
is classically forbidden by the energy barrier. This is 
very general, even for multidimensional non-integrable 
systems: the splitting AE n = \E~ — E+\ between some 
nearly degenerated doublets provides a quantitative ma- 
nifestation of the tunnelling between the phase-space re- 
gions where the corresponding eigenstates are localised. 

Rather than computing the poles of E h-> G(a, —a, E), 
another systematic strategy to obtain one individual 
splitting [48| is to start with Herring's formula [25l - [27j 
that relies on the knowledge of the eigenfunctions outside 
the classically allowed regions in pha se-space. Here, we 
will propose alternative formulae [28[ that involve traces 
of a product of operators, among them the evolution ope- 
rator, 

U(T) d ^ e- iiiT / h 

OO 

= £ 0»> <#l ^ EtT ' h + l*n) (4>n\ e- iE " T/ft ) , (5) 
n=0 

analytically continued in some sector of the complex 
time domain. This approach, which privileges the time- 
domain, provides, of course, a natural starting point for a 
semiclassical analysis in terms of classical complex orbits. 

The simplest situation occurs when the tunnelling dou- 
blet is made of the two lowest energies Eq < Eq = 
Eq + AEq of the spectrum. If hco denotes the energy dif- 
ference between the nearest excited states and Eq (for 
the double well potential Ef ~ Eq + fku where ui is the 
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classical frequency around the stable equilibrium points) , 
by giving a sufficiently large imaginary part to — T, 



-cjIm(T) > 1 , 



(6) 



we can safely retain the n = terms only, which exponen- 
tially dominate the trace of d5j . To be valid, this approx- 
imation requires that we remain away from a quantum 
resonance where the definition of the tunnelling doublet 
is made ambiguous by the presence of a third energy level 
in the neighborhood of Eq and Eq . Then we have im- 
mediately 



tan 



( TAE \ ti(SU(T)) 
V 2ft J ~ tr(U(T)) 



(7) 



When the tunnelling splitting is smaller than fko by sev- 
eral orders of magnitude, we can work with a complex 
time such that 



\T\AE, 



21) 



< 1 



(8) 



remains compatible with condition © and, therefore, 



AE Q ~ A (T) ^ 



s 2h tr(SU(T)) 
' W t r(£7(T)) 



(9) 



will provide a good approximation of the tunnelling split- 
ting for a wide range of T. Though condition (j8]) is widely 
fulfilled in many situations this is not an essential condi- 
tion since one can keep working with tan" 1 . Numerically, 
the estimation ^ has also the advantage of obviating a 
diagonalization. 

When we want to compute a splitting, due to tun- 
nelling, between an arbitrary doublet, the selection of 
the corresponding terms in the right hand side of ([5]) 
can be made with an operator TI n that will mimic the 
(a priori unknown) projector + \(j>~) It 

will be chosen such that its matrix elements are localised 
in the regions of phase-space where 4>^ (p, q) are domi- 
nant. Under the soft condition \T\AE n /(2h) <C 1, we 
will therefore take 



AE n ^A n (T)^^^$4^ 



IT tr(n n U(T)) 



(10) 



The localization condition on the matrix elements of tl n 
is a selection tool that replaces ©. In that case, there is 
a battle of exponentials between the exponentially small 
matrix elements (</>™| n n and the time dependent 
terms e~^ Em ~ En ^ T ' h . The last term would eventually 
dominate for the lower energy states (E m < E n ) if — Im T 
could be increased arbitrarily. But once a ImT < is 
given, one expects to recover a good approximation of the 
excited splitting by increasing Re T since, when Re T ^> 
| ImT|, we actually recover the real time case. 

The next step consists in computing A n {T) by semi- 
classical techniques and, then, we will add some more 



specific prescriptions on the choice of T, in order to im- 
prove the accuracy of AE n . We will illustrate how this 
works in the examples of sees. |lVl |V] and IVII 

Let us mention another way to select excited states, 
that we did not exploit further. With the help of a posi- 
tive smooth function F(u) that has a deep, isolated min- 



imum at u — 0, say F(u) 



,2N 



with N strictly positive, 



we can freeze the dynamics around any energy E by con- 
sidering a new Hamiltonian H'(p,q) =F(H(p,q) - E)). 
Classically, the phase-space portrait is the same as the 
original one obtained with Hamiltonian H except that 
the set of points H (p, q) = E now consists of equilibrium 
points. The quantum Hamiltonian H' =F(H — E) has 
the same eigenfunctions as H but the corresponding spec- 
trum is now F(E^ — E). By choosing E = £"+, the 
doublet i?„ yields to the ground-state doublet AE' n — 
F(AE n ) and we can use the approximation §§§ with U' = 
cxp(-iH'T). With the F given above, we have 



AE„ 



[A;(T)] i/(2A0 



def 



2fttr(S U'(T)) 
iT tr (U'(T)) 



1/(2N) 



(11) 



III. SEMICLASSICAL EXPRESSIONS 

A. Hamiltonian dynamics with complex time 

Formally, the numerator and the denominator of the 
right hand side of (|9]) can be written as a phase-space 
path integral of the form 



e iS[p, q ;t]/h D [p] D [ g ] 



(12) 



where the continuous action is the functional 



S[p,q;t] d = (p( s )j^ 



dt 



(s)-H(p(s),q(s))-(s)jds. 

(13) 

The subset V of phase-space paths, the measure D[p]D[q] 
and the action (fl3|) appear as a continuous limit (r — > 0) 
of a discretized expression whose precise definition de- 
pends on the choice of the basis for computing the traces 
but, in any cases, involves a typical, finite, complex time 
step r (see also appendix [X| . The complex continuous 
time path s t(s) is given with fixed ends f (sj) = tj = 0, 
t(sf) = tf = T. Because the slicing of T in small complex 
time steps of modulus of order r is arbitrary, the inte- 
grals of the form (|T2|) remain independent of the choice 
of t(s) for Si < s < sf as long as Imt(s) is non-increasing 
in order to keep the evolution operators well-defined for 
any slice of time [l9j]. In the following we will denote 
by [t] such an admissible time-path. In a scmiclassical 
limit (keeping the order lim^o limT^o 0)i the dom- 
inant contributions to integrals (|12p come from some 
paths in V that extremise S, i.e. from some solutions 
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of Hamilton's equations 



known Van Vleck approximation for the propagator: 



dp 
ds" 
dq 
ds 



dH dt 
dq ds ' 
dH dt 
dp ds 



(14a) 
(14b) 



with appropriate boundary conditions imposed on some 
canonical variables at Sj and/or Sf. When H is an an- 
alytic function of the phase-space coordinates (p,q), we 
can take the real and imaginary parts of equations (1141) . 
use the Cauchy-Riemann equations that render explicit 
the entanglement between the real part and the imagi- 
nary part of any analytic function f(z) : Re(d//dz) = 
<9(Re/)/<9(Rez) = 9(Im f)/d(lm z) and Im(d//dz) = 
d(lmf)/d(Rez) = —d(Ref)/d(Imz), and then obtain 
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Re 
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d Reg 


d(— Imp) 
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Re 
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dim q 


dReg 
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Re 


ds 
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dlmq 


d 




ds 


d{— Imp) 



r di 

ds 
r dt 
ds 
r dt 
ds 

Re ( H 



.dt 



(15a) 
(15b) 
(15c) 
(15d) 



Therefore, the dynamics described in terms of complex 
canonical variables is equivalent to a dynamics that re- 
mains Hamiltonian — though not autonomous with re- 
spect to the parametrisation s whenever dt/ds varies 
with s — , involving twice as many degrees of freedom 
as the original system, namely (Re q, Im q) and their re- 
spectively conjugated momenta (Rep, — Imp). The new 
Hamiltonian function is then Ke(Hdt/ds) but it would 
have been equivalent (though not canonically equivalent) 
to choose the other constant of motion Im(Hdt/ds) as a 
Hamiltonian. 

What will be important in what follows is that [t] 
will not be given a priori. Unlike in the standard in- 
stanton approach where t is forced to remain on the 
imaginary axis, we will see that for describing tunnelling 
it is a more efficient strategy to look for some com- 
plex paths (p(s),q(s)) that naturally connect two phase- 
space regions and then deduce [t] from one of the equa- 
tions (TT4")) . It happens that in the two usual textbook 
examples that we mentioned in the first paragraph of 
our Introduction, the complex time [t] has a vanishing 
real part, but in more general cases where tunnelling be- 
tween excited states is studied, this is no longer true. 



B. Trace formulae 

Let us privilege the ^-representation and consider the 
analytic continuation for complex time T of the well- 



G(q f , qi ;T) ~ V(-ir° 
h—tO * — ' 



I det 



d 2 S e 



2irh dqidqf 



jS /h 



(16) 

The sum involves (complex) classical trajectories o in 
phase-space i.e. solutions of equations ([T5|) with q(si) — 
Qi, <l{ s f) — Is f° r a gi ven [t] such that f(sj) = 0, 
t(sf) = T. The action S is computed along o with defi- 
nition (fT5T) and is considered as a function of (qf,qf,T). 
The integer k encapsulates the choice of the Riemann 
sheet where the square root is computed; it keeps a record 
of the number of points on o where the semiclassical ap- 
proximation (|16|) fails. As far as we do not cross a bifurca- 
tion of classical trajectories when smoothly deforming [t], 
the number of o's, the value of S and n do not depend 
of the choice of [t]. The numerator (resp. denomina- 
tor) of Ao(T) are given by the integral J G(riq,q;T)dq 
where 77 = —1 (resp. r\ = +1). Within the semiclassi- 
cal approximation, when contributes to the propaga- 
tor G(?7g,q;T), we have to evaluate 



(-!)"•/ dg 



det 



d*S n 



2i:h dqidqf 



JSo(vq,q,T)/h 



(17) 

The steepest descent method requires the determination 
of the critical points of S (r)q,q,T). Since the momenta 
at the end points of are given by 

Pi = p(si) = -d qt S (qf,qi;T) ; (18a) 
Pf = P( s f) = 9 qf S„(q f ,qi;T) , (18b) 

the dominant contributions to tr([/(T)) come from perio- 
dic orbits, i.e. when (pf,qf) = (jPi,q%), whereas the dom- 
inant contributions to ti(SU(T)^j come from half symme- 
tric periodic orbits, i.e. when (pf,qf) = {—pi,—qi). As 
explained in detail in the appendix \^ one must distin- 
guish the contributions of the zero length orbits e (the 
equilibrium points) from the non-zero length periodic or- 
bits 0. For one degree of freedom, their respective con- 
tributions are given, up to sign, by 



-xH(jp e ,q e )T/H 



3 A c T/2 



rje~ 



-A c T/2 



(19a) 



with ±A C the two Lyapunov exponents of the equilibrium 
point e = (p e ,q e ), and 



(-1) A 



\/—2rjmh 



dE 
dT 



jSJh 



(19b) 



where is a periodic orbit of period T = T if rj = +1 
and half a symmetric periodic orbit of half period T = T 
if r) = — 1. The sum runs over all the branches j3 that 
compose the geometrical set of points belonging to 0. Tp 
is the characteristic time (|A12[) on the branch (3. The 
energy E is implicitly defined by (|A10[) and /i essen- 
tially counts the number of turning points on 0. 
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Expressions (p~9|) are purely geometric; their classical 
ingredients do not rely on a specific choice of canonical 
coordinates and they are independent of the choice of the 
basis to evaluate the traces. 

In the general case, the most difficult part consists in 
determining which periodic orbits contribute to the semi- 
classical approximation of the traces. Condition T = T 
is necessary but far from being sufficient; the structure 
of the complex paths (keeping a real time) may appear 
to be very subtle [29| and have been the subject of many 
recent delicate works [30[ • Even for a simple oscillating 
integral, the determination of the complex critical points 
of the phase that do contribute is a highly non-trivial 
problem because it requires a global analysis: one must 
know how to deform the whole initial contour of integra- 
tion to reach a steepest descent paths. 

Our strategy consists in retaining the terms (|19b[) for 
which we can choose the complex time [t] to constrain 
the (half) periodic orbit o to keep one of the canonical 
coordinate (say q) real. This o surely contributes because 
we do not have to deform the q part of the integration do- 
main of (1121) in the complex plane. But we will see in the 
next section that two different periodic orbits with real q 
correspond to two different choices of [t] : For a chosen [t] , 
only some isolated points, if any, in phase space, will pro- 
vide a starting point (j>i, qi) of a trajectory with q(s) real 
all along o. When sliding slightly the initial real coor- 
dinate qi in the integral ([T7|) . it requires to change [t] 
as well to maintain q(s) real on the whole o. This is 
not a problem since all the quantities involved in expres- 
sion (fT7|) are [t] -independent if the deformation of [t] is 
small enough not to provoke a bifurcation of o, that is, 
whenever the initial point does dot cross one of the turn- 
ing points, which are the boundaries of the branches (3. 
The integral (|I7[) may be computed with a fixed time- 
path or an adaptive one for each separated branch. This 
computation, presented in the appendix [XI generates the 
sum over all branches that appears in formula (fl9b|) . To 
sum up, to keep a contribution to the traces, we must 
know if we can choose a shape of [t] in order to pick a 
real-g periodic orbit with a specific qi. 

This construction is certainly not unique (one may 
choose other constraints) so we will use the intuitive prin- 
ciple that the periodic orbits we choose will connect the 
two regions of phase-space that are concerned by tun- 
nelling i.e. where 4>n are dominant; this is justified by 
the presence of the operator TI n in the right hand side 
of (JIOJ). The prcfactor n n (q,q') d = (q\ 77„ \q') will select 
in the integral 

J ' &q&q l n n {q,q')G{ m l ,q:T) , (20) 

a domain around the projection onto the g-space of the 
appropriate tori. In the specific case of the ground-state 
splitting, condition ^ does the job of TIq: the large 
value of ImT requires that the orbit approach at least 
one equilibrium point and then follow a separatrix line. 



In general we will not be able to prove that other com- 
plex periodic orbits give sub-dominant contributions but 
the examples given in the following sections are rather 
convincing. Moreover, our criterion of selection allows 
us to justify the four rules presented in [3f|, sec. IIA] 
for computing the contributions of complex orbits to the 
semiclassical expansions of the energy Green's function. 

IV. APPLICATION TO THE DOUBLE- WELL 
POTENTIAL 

Let us show first that our strategy leads to the usual in- 
stanton results for a Hamiltonian of the form ((IJ with an 
even V having two stable symmetric equilibrium points 
at q = ±a (Fig. []Ji); in their neighbourhood, the fre- 
quency of the small vibrations is uj. We are interested 
in the ground-state splitting AEq for h small enough in 
order to have AEq <C fow. 

Before we make any semiclassical approximations, we 
check in Fig. [2] the validity of the estimation © on the 
quartic potential by a) verifying that Aq(T) is almost 
real and independent of the choice of the complex T pro- 
vided that conditions © and © are fulfilled, and b) by 
checking that this constant gives a good approximation 
of the "exact" AEq computed by direct numerical diag- 
onalization of the Hamiltonian [32[. 

As explained at the end of section sec. IIIIl in phase- 
space we will try to find some time-path [t] that allows 
the existence a (half) symmetric periodic orbit o that 
connects two tori at (real) energy E > in the neigh- 
bourhood of (p, q) = (0, ±a) while q remains real. If we 
impose t(si) = and t(s/) = T, then E = E (qf,qi,T) 
is implicitly given by relation (|AI 0[) . We will denote 
by q r {E) (resp. q' r (E)) the position of the turning 
point at energy E > that lies in between q = 
and q — a (resp. that is larger that a). The two 
branches p± (g, E) = ± y/2{E — V(q)) are either purely 
real when V(q) E or purely imaginary when V(q) ^ E. 
Then, from equation (|14b[) . 



is purely real in the classically allowed region, while 
purely imaginary in the forbidden region. Therefore the 
complex time path [t] must have the shape of a descend- 
ing staircase whose steps are made of pure real or pure 
imaginary variations of time (see Fig. [3]). 

The complex orbit with real q can be represented in 
phase-space as a continuous concatenation of paths, fol- 
lowing the lines E — (Rep) 2 /2 + V[q) in the allowed 
region and — E = (Imp) 2 /2 — V(q) in the forbidden re- 
gion. It is natural to represent o in the three dimensional 
section Iraq = of the complex phase-space with axes 
given by (Re q, Rep, Imp): the junctions at the turning 
points lie necessarily on the (Rep = 0, Imp = 0) axis (see 
Figs-Hfc and d). A periodic orbit is made of a succession 
of repetitions of 
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Figure 1: (colour online) In the case of a double well potential 
shown in a), tunnelling can be described by a solution of (|14[) 
with a purely real q evolving from one well to the other. In 
phase-space, this trajectory appears to be a concatenation of 
two types of curves that join on the (Reg)-axis: 1) a trajec- 
tory that lies in the phase-space plane (Reg, Rep) with real 
variations of t at energy E with the potential V and 2) a 
trajectory that lies in the phase-space plane (Reg, Imp) with 
imaginary variations of t at energy —E with the potential —V 
shown in b). In c) a family of constant energy curves is shown 
(horizontal green (light gray) for the first type, vertical red 
(dark gray) for the second type). In d) for a given energy, we 
show how the three curves [, c, r glue together at the turnin g 
points (0, ±g r ). & denotes the separatrix Imp = ±y/2V(q). 



(i) primitive real periodic orbits r with energy E in the 
right region, that is, such that q r {E) < q{s) < q' r (E) and 
Imp(s) = 0; 

(ii) primitive complex periodic orbits c in the central re- 
gion with purely imaginary p and real g, that is, such 
that —q r (E) < q(s) < q r {E) and Rep(s) = 0; 

(iii) primitive real periodic orbits [ with energy E in the 
left region. They are obtained from the periodic orbits r 
by the symmetry S. 

By denoting T t (E) and T C (E) the (real, positive) peri- 




b) 

0.25, 



0- 



-0.25 
-0.5 
-0.75 
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Figure 2: (colour online) For V(q) = (q 2 — a 2 ) 2 with a — 1 
and ft ~ 1/12, we have plotted in the upper graph a) the real 
part of Ao defined by ([9| as a function of the complex T. 
It becomes constant for ImT large enough for having ([6]) 
(1/w = l/(2V2o) ~ .35) regardless of ReT. In the same 
range, Im (A () (T)) is negligible compared to AEo ~ 4.4 1CT 10 
and condition (8) is largely fulfilled. In the lower graph b), we 
have computed the "exact" value AEo by direct diagonaliza- 
tion. The black thick solid line corresponds to ln(A£o) from 
which we have substracted A (ft) = f — S c (ftw/2)/(2ft)+ln(ftw/7r) 
in order to emphasize the contribution of the prefactor. The 
blue solid line corresponds to In | Ao| — A(ft) for T ~ — 4i and 
become indistinguishable from the previous one for 1/ft > 3. 
With the same substraction, our first proposal (|39[) gives the 
constant (In y/n) ~ 0.6 and our second one (|44L which is 
the same as the formula in Landau and Lifshitz [331 . sec. 50, 
problem 3], corresponds to (dotted line). Garg's formula [27l . 
eq.1.1] gives the constant (In \/ir/e ~ 0.07) (dash-dotted line). 

ods of the primitive periodic orbits r and c at energy E 
respectively, we have 

T {E) = w,T z {E) - iw c T c {E) - T C (E) (22) 

with the winding numbers w t and w c being non-negative 
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integers. The trajectory o may contribute to the deno- 
minator (j] = +1) or to the numerator (77 = —1) of the 
right hand side of ((9|) provided that T — T. For 77 = +1, 
is periodic whereas, for 77 = —1, is half a symmetric 
periodic orbit. We will keep the contributions of all or- 
bits for which a staircase [t] can be constructed. We can 
understand from figure [3] that orbits differing by a small 
sliding of their initial qi can be obtained by a small mod- 
ification of the length or height of the first step. All these 
contributions are summed up when performing the inte- 
gral (|17[) on one branch f3 and correspond to one term 



in the sum Y2p m formula (119b[) . When at least one 
of the winding number is strictly larger than one, several 
staircase time-paths can be constructed while keeping re- 
lation (f!2!2)) : they differ one from each other by a different 
partition into steps of the length and/or heights of the 
staircase. The corresponding orbits can be obtained 
one from each other by a continuous smooth sliding of 
the steps of the staircase but during this process one 
cannot avoid an orbit starting at a turning point where 
a bifurcation occurs. 




T c /2 



T% — T 



T 




Figure 3: (colour online) For the double well in Fig. [T] at a given energy < E < Vmax, when T = 2T t (E) — iT c (E)/2, for each 
one-step time path, there exists a unique half symmetric periodic orbit with energy E keeping Iraq = starting from the left 
allowed region. Its initial point is uniquely defined: For the time-path depicted in the left column, (pi, qi) is such that it takes 
exactly the real time r to reach the turning point (0, — q r ) following the primitive orbit [. Then it follows half of the orbit c and 
ends in the right well winding once along r (a). When starting at q- < qi, there still exists a half periodic orbit joining (j4, q-) 
to (—p'i,—qi) in time T for some p\ but unlike the previous one it gets outside the three dimensional section Img = (its 
projection is schematically represented by the dashed line in a'). However an appropriate shift of the time step to r' > r 
(central column) turns the latter smoothly into a trajectory with Imq = (b'). During this process, the first trajectory starting 
at qi cannot keep a purely real q anymore (b). If we shift the step by exactly T r (right column), we recover a trajectory starting 
at qi with purely real q (c) but with a different topology since it is now winding once along [. The change of topology occurs 
when the initial position reaches a turning point while sliding the step. 



In the right hand side of (fT6)) , the sum involves sev- 
eral trajectories differing one from each other by the se- 
quence of turning points that are successively encoun- 
tered along 0. Therefore, to compute the dominant con- 
tribution of the non-zero length orbits to the numerator 
and to the denominator of (j9]) , we will add all the contri- 
butions of the topological classes of orbits, each of them 
uniquely characterised by an ordered sequence of turning 
points [pi , P2 , . . . ] , in other words by a partition of w v 



and w c into integers and by the branch /3 where its star- 
ting point lies. We can therefore express our result in a 
way that can be applied to cases more general than the 
double well: the total contribution of the non-zero length 
paths to the numerator (77 = —I) and to the denominator 
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(r? = +1) of © is 

« [P1.P2,...] V ' 

s denotes a section of an energy surface in complex phase 
space corresponding to one purely real canonical variable. 
\pi, p2, ■ ■ . ] is an ordered sequence of (not necessarily dis- 
tinct) turning points that belong to a section s. The sum 
concerns all s (different energies may be possible) and 
[pi, p2, ■ ■ . ] such that we can construct on s, with an ap- 
propriate choice of [t], a periodic orbit o if r) = +1 (a half 
symmetric periodic orbit o if r; = —1) of period T = T. 
The branch j3 is the one where o starts, the sequence of 
turning points that are successively crossed by o is ex- 
actly [pi,p2, ...]• 

In the case of the double- well, for an energy be- 
low K nax , a section s for real q has four turning 
points (0, ±g r ) and (0,±q' r ). Only the points on s such 
that — q' r < q < q' r can provide starting points of a perio- 
dic orbit. They belong to one of the three closed loops r,c, 
[ that connect on the axis (p = 0, Imq = 0) at the turning 
points (0, ±<7 r ). Once T is given, the condition T = T 
will select a finite set of energies for which we have (|2"2")l 
with integers w t and w c . Any o with winding numbers w t 
and w c will have an action given by 

S = w x S z + iw c S c + i 1 - - S c (24) 

and an index 

1 77 

Mo = w v + w c H — (25) 

(in the real case, 2/x„ is computed in the same way as the 
Maslov index: it counts the number of turning points 
encountered along 0). These quantities are independent 
on the choice of the six possible starting branch (3 (each r, 
c, [ is made of two branches). When the orbit starts 
on r or on I, we will have Tp — T x /2 and when starts 
on c, we have Tp = — iT c /2. 

For — Im T larger than the oscillation period in the 
central well of — V, w c — 1 is the minimum value of the 
winding number when 77 = +1 (just one back and forth 
trip around c) while it is w c = when rj = — 1 (just one 
half of c is concerned). For these orbits, condition © 
forces c to stay near the separatrix 6 defined by Imp = 
±y/2V(q) and, thus, r must lie in the immediate neigh- 
bourhood of the equilibrium point (p,q) = (0,a). These 
orbits will give the dominant contribution because they 
have the smallest Im S among all the other possible or- 
bits involving repetitions of c. Indeed, for a fixed T, all 
the orbits that may contribute semiclassically are such 
that ImS = -2Im(T)S c /T c and S c /T c is a decreasing 
function of the energy E = — E when c is inside the 



separatrix since 

d (S C (E)\ = 
dE \T C (E)J 

4 dT (E) r qr( - E) 

Therefore Im S reaches its minimum when E is at 
its maximum, that is E — > + . The only possible 
equilibrium point contributing to tr(St/(T)) is the ori- 
gin (p/,qf) = {~Pi,qi) = (0,0); it is also subdominant 
because H(0, 0) = V(0) = Knax (which can be seen as 
the limit of S C (E)/T C (E) when E -> -F max ) is strictly 
larger than S C (E)/T C {E) for E > -V max . 

Assuming that only orbits with real q do contribute 
to tr(S U(T)), we have proven that the dominant contri- 
bution is given by the half symmetric orbits at energy 
E such that T C (E) = -2ImT (iu c = 0). However, for 
such an orbit to exist we cannot choose ReT arbitrarily 
since it must be an integer multiple of T t (E). To put 
it in another way, for a given E, we will choose T = 
w t T t (E) — iT c (E)/2 such that an orbit with a real 
q exists. Condition © will be fulfilled if E is suffi- 
ciently small. We must now enumerate all the topological 
classes concerned by the sum (f23|) : When starts on the 
upper branch Imp > of c, it will reach the turning 
point p\ = (0, q r ) then wind times around r alterna- 
tively crossing (0,q' r ) and (0,q r ) before turning back on 
the lower branch Imp < of c. Starting on the lower 
branch of c corresponds to the symmetric trajectory and 
provide the same contribution with Tp = —iT c /2. When 
starting on the upper branch of r, crosses p\ = {0,q' r ) 
first, then reaches P2 = (0, q r ). Then it may go on wind- 
ing r times around r then take the lower branch of c up to 
the turning point (0, — q r ), wind w t — 1 — r times around [ 
and eventually join the symmetric of its starting point on 
the lower branch of t. There are exactly w x such topo- 
logical classes because we can take r = 0, w t — 1. If 
we start from the lower branch of r or on one of the two 
branches of [, we obtain the same contribution and ex- 
haust the possible topological classes. The sum of Tp on 
all P's and classes is then 2(-iT c /2) + 4w r T r /2 = 2T and 
keeping only the w c = solutions, the sum (|23p reduces 
to 

with being one half symmetric orbit of energy E de- 
fined implicitly by 

T = w z T x {E) - iT c (E)/2 (28) 
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with 



T t (E) 



T C (E)=4 



dq 



lAE) ^2(E-V{q)) 
dq 



q r (E) 



l2(V(q) - E) 
We have S (E) = w z S x {E) + iS c (E)/2 with 



(29a) 
(29b) 



S t (E) + ET t (E) = S t (E) = 2 / J2(E - V(q)) dq, 

lq T {E) V 

(30a) 

r q r (E) 

2(V(q) - E) dq. 

(30b) 



q'r(E) 



clef 



S C {E) - ET C (E) = S C (E) = 4 



The dominant contributions to tr([/(T)) comes from the 
two stable equilibrium points e = (0, ±a) for which A c = 
iui. The contribution of (0, 0) is sub-dominant as well as 
the contribution of any periodic orbit o which necessarily 
turns around about c during Re T then follow an orbit c 
near 6 during — Im T before coming back to its initial 
point. Together with ©, 77 = 1, the two stable equili- 
brium points give two identical contributions (I19a[) and 
we have 



tr(U(T)) ~ 2e 



-iuiT/2 



(31) 



which of course could have been deduced directly 
from Eq ~ hui/2. Collecting all these results in the right 
hand side of ([5]), we obtain 



2ift f dT, 



A ° hZo V 7T [ dE 



-1/2 



x exp (-±_(§ C (E) + (E- tuv/2)T c (E)^j 



' ( ' X ' > ' ~7T~ v — (E — hw/2)T z {E) — frnj J . (32) 

When E —> + , we have the following asymptotic ex- 
pansions (see appendix FBI) 

(33a) 
(33b) 

(34) 
(35) 



S t (E) 
with 



2-kE 



^def 



BE 2 + o(E 2 ) ; 



FA- 



1 



dq ; 



B 6* n 



\y/2V(q) a-q f 



(the superscript in parenthesis indicates the order of the 
derivative of V). The differentiation of expressions (|3"3")l 
with respect to E leads to the asymptotic expansions 
for — T c and T r . From the first one we can extract the 
exponential sensitivity of E on ImT: 



E~2aVe 2 V IraT 



(36) 



From relation (|2"8j) we can see that 2dT /dE = 
2w t dT t /dE — idT c /di? is dominated by the last term 
if E is small. 



dT 
dE 



1 dT c 

2 ~dE 



1 

ujE 



(37) 



Inserting all these asymptotic expansions in the right 
hand side of (|3"2")l , we get 



A ^ J 4ha2 U S c A c E/(tm) c -g t "n .2/,. 



fi-vO V 7T 



x e 



-\w t BE(E—huj) /h 



(38) 



This expression can be turned into the usual jwkb ex- 
pansion exp (a (E)/h + ai(E) In ft + a 2 (£') + o(l)): As 
soon as we have condition ([6]), from expression (|36j) we 
see that E is exponentially small and we obtain a (0) = 
— S c (0)/2 and ai(0) = 1/2. To obtain the correct value 
of et 2 , we must proceed to fine-tune the choice of T. A 
criterion is to impose on Ao to have a vanishing imag- 
inary part at any order in ft consistent with the jwkb 
expansions used so far. From (|32[) . we will choose E 
such that S t (E) = fm, which is exactly the Einstein- 
Brillouin-Keller quantization condition for the ground 
state in one well. This leads to E = haj/2 + o(ft). 
Then T c = -2 In (ft/(4a 2 u;)) /lj + AA/uj + o(l) and 



A 



-S c (hu/2)/(2h) 



2auj 



ehuj 



e A e -s c (0)/(2h) 



_ (39) 
which differs from [23, eqs. (1.1) and (1.4)] by a reason- 
able factor -y/e ~ 1.6. This discrepancy, already noticed 
in [3J, sec. V], which appears in the third order term in 
the ft-expansion, comes from the different kind of approx- 
imations involved in our approach on the one hand and 
in Herring's formula on the other hand. 

We are also able to obtain a formula for the splitting 
of the excited states that is consistent with the result 
given in [27], eq. (Bl)]. Using a semiclassical approxima- 
tion for the matrix element of 7T„, we explain in detail in 
appendix [Cl how to obtain A n (T). For one dimensional 
systems whose energy surface E n is made of two branches 
(two Riemann sheets in the complex plane), we can in- 
sert (|C8p into (1C6[) and get one of the main result of this 
paper: 



A„(T) 



h-^o 2T ^ 

[P\,P2,-. 



(_ iY° +1 e^°{ b fi'> b P, E -»-)/K . (40) 



To see how formula f|40[) works in the case of the double- 
well potential, we choose the quasi- mode localised 
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on the right torus r at energy E n . This torus is made of 
two branches labelled by the sign of p and we can choose 
a common base point for these two branches, namely 
b± = q r (E). On the symmetric torus [, the two base 
points will be b± = —q r (E). Then S (bp> ,bp , E n ) and 
the index /j, do not depend on the choice of the initial 
and final branch. The orbits o that go from r to [ must 
correspond to a T of the form 

1 



T = t + w^T x (E n ) - i [w c 



- ) T c (E n ) (41) 



for non- negative integer (w c , w t ) and a fraction of time r 
strictly smaller that T c (E n ) that depends on the initial 
and final conditions (those are not necessarily symme- 
tric). Then we have 



S o (b^,b ,E n ) = w v S v (E n )+i[ w c + -) S c (E n ) (42) 



and we take 



W c 



1 



(43) 



For the same reason as previously explained the domi- 
nant contributions come from those orbits where w c = 0. 
In order to mimic a real T, we will choose large wind- 
ing numbers w z such that ReT = w t T t (E n ) ^> ImT = 
T c (E n )/2. Because of the quantization condition in the 
right well S z (E n ) = (n + 1/2)2ttH, the rapid oscillations 
exp(iRe S /K) disappear (or inversely if we want to main- 
tain A„ real to first order, we recover the usual quanti- 
zation condition). Then we obtain 



21i 



-S c {E n )/(2h) 



H^o T,{E n ) 



(44) 



The classical frequency l/T t (E n ) attached to r is es- 
sentially of order oj/(2it). A„(T) becomes independent 
of T for large w z : the behaviour of T is mainly gov- 
erned by w t T t , then the 1/T prefactor in (14"0|) is compen- 
sated by the increasing number of identical terms in the 
sum since we have seen that the number of topological 
classes of orbits increases linearly with Aw t (the factor 4 
comes from the two possible initial branches /3 = ± and 
the two possible final branches /3' = ±; in other words 
from the sequences of [p\, p 2 , . . .] beginning either by q r 
or q' r and ending either by — q r or —q' r ). The discrep- 
ancy between (|44l) and Garg's formula is just the fac- 
tor g n given by [27], eq. (B2)] (see also[3^ eq. (3.41)]) 
that tends to 1 when n increases: go — x/n/e ~ 1.075, 
gi ~ 1.028, . . .There is also a ratio of order one, more 
precisely \f2pK ~ 0.8, between estimations (|39|) and (|44"]) 
taken for n = 0; the second is slightly better and coin- 
cides with the formula given in Landau and Lifshitz (33l . 
sec. 50, problem 3]. Again, these discrepancies come from 
the different nature of the approximations that arc in- 
volved. 

Let us end this section by a short comment on the con- 
nection with the usual instanton theory where Re T = 



and where ImT — > — oo. This regime, that allows to 
select the ground-state doublet only, is included in our 
approach because the instanton trajectories appear to be 
the limit of c getting closer to the separatrix © whereas 
the classical real oscillations in the wells shrink to the 
equilibrium points. All along this paper we emphasize 
that the phase-space representation is particularly ap- 
propriate and it is straightforward to recover the usual 
picture of instantons (for instance q(jt) versus it) from 
our figure [T]d). 



V. DYNAMICAL TUNNELLING FOR THE 
SIMPLE PENDULUM 

The simple pendulum corresponds to V(q) = — 7COsg 
with 7 > and strictly periodic boundary conditions that 
identify q = — tt and q — ir. At energy E > 7, the classi- 
cal rotation with p = \J2(E — V{q)) can never switch to 
the inverse rotation with — p = —y / 2(E — V(q)). At the 
quantum level, the Schrddinger's equation for the sta- 
tionary wave-function (j) 



7Cosg 



0(g) = E<t>{q) 



(45) 



leads to the Mathieu equation [36j |. 

y"(x) + (a - 2gcos(2x))y(x) = (46) 

with x = q/2, y(x) = 4>{2x), a = 8E/H 2 and g = 
—A'j/h 2 . The 27r-periodicity of 4> forces y to be ir- 
periodic. The eigenfunctions can be classified accord- 
ing to the parity operator: Even 7r-periodic solutions 
exist only for a countably infinite set of characteristic 
values of a denoted by {d2 n } with n = 0, 1, . . . Odd solu- 
tions correspond to another set, {62^}, with n = 1, 2 . . . 
(the {a2«+i, &2n+i} correspond to 7r-antiperiodic solu- 
tions and will be rejected). The discrete energy spec- 
trum corresponding to even and odd solutions of (|45[) 
is then {E+ = h 2 a 2n /8,n = 0,1...} and {E~ = 
h 2 b2n/8,n = 1,2...} respectively. Any eigenstate |</>^) 
with energy E^ ~ E n > 7 has its Husimi distribution 
spread symmetrically between the two half phase-space of 
positive and negative p, near the lines ± y / 2(E n — V(q)) 
that define two disconnected tori in the cylindrical phase- 
space. If we prepare a wave-packet localised on the line 
p = yj2(E n — V(q)), since it is no longer a stationary 
state, its average momentum will oscillate between two 
opposite values, with a tunnelling frequency equal to 
AE n /h where 



El 



e: 



h 2 

-c7\ a 2n 



b 2 n) 



(47) 



is the splitting between the two quasi-degenerate eigenen- 
ergies. To compute AE n , we will use ([TU]) with the ope- 
rator II n very much like the exact projector |<^+) (</>+| + 
\4>n) {4>nV ^s matrix element (p'\ II n \p) will vanish 
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rapidly as soon as p or p' lie outside the region of the 
two tori. The main contribution to the semiclassical ex- 
pansion of tr(S IJ n U(T)) will come from classical trajec- 
tories that connect two symmetric tori. The trace will 
be semiclassically computed in the momentum basis and 
we will choose the complex time path to maintain p real. 
To construct one trajectory at energy E > 7 connecting 
the two tori requires to have purely imaginary q when- 
ever — p r < p < p r with p r = y/2(E — 7) being the clas- 
sical turning point in momentum. More precisely, q is 
given by cosg = — 1 — (j% — p 2 )/(2j). From (|14a[) 



d* 
ds 



1 dp 

7 sin q ds 



(48) 



we see immediately that sin (g(s)) and dt/ds will be 
real when p(s) > p r or p(s) < —p T and purely 
imaginary otherwise. In the latter case, since we 
keep P2(s) =lmp(s) = 0, equations (|T5b .c) lead to two 
possible families of solutions (i) qi(s) = Req(s) = and 
(ii) <7i(s) = ±7r. Then, with pi = Rep and q2 = — Img, 
equations (fl"5"a.d) become 



dpi , . (At 



dq 2 



= Pi 



dt\ 



(49a) 
(49b) 



and are associated with a real time dynamics governed 
by the Hamiltonian 



H(pi,q 2 ) =Pi/2=F7chg 2 



(50) 



with "— " corresponding to case (i) and "+" correspond- 
ing to case (ii). Then, instantons correspond to trajecto- 
ries evolving in the transformed potential V(q) = ^7 chg 
rather than the usual inverted potential — V(q). We 
will choose a one step complex time path as in Fig. [3] 
and we will represent the orbits in a three dimensional 
space (Re g, Rep, Im q) and the connection between the 
allowed and the forbidden trajectories occurs on the 
axes (Reg = 0,lmg = 0) and (Reg = ±7r,Img = 0) 
(Fig. 2]a,). Trajectories in family (i) escape from the un- 
stable point at Rep = 0,lmg — without coming back 
to the plane Rep = 0,Reg = 0. Only orbits in family 
(ii) can be used to produce periodic orbits. A typical pe- 
riodic orbit connecting two symmetric rotations of the 
pendulum is given in Figj4]3). We will choose ImT to be 
precisely the half period of the periodic orbit c of fam- 
ily (ii) at energy E n and Re T to be (an integer multiple 
of) the period of rotation of the pendulum at the same 
energy. This choice exhibits the dominant contribution 
to tr(S fl n U(T)j . We can reproduce the same reasoning 
that led to (|44p with the role of p and q being exchanged. 
Now T r = dS z /dE is the typical frequency on the torus 
at energy E n . 

Keeping only the contributions that provide a positive 




Figure 4: (colour online) With the choice of a complex- 
time path given in [3| (one-step path), the contributions 
to tr(S/7„ (7(T)) come from one half symmetric periodic or- 
bits (b) that are a composite of 1) (repetitions of) a periodic 
orbit r of period T c in the phase-space plane (Reg, Rep) at 
energy E > 7 (recall that the planes q = ±7r are identified); 
2) one half of a periodic orbit c of period T c in the phase-space 
plane (Img, Rep) at energy E. 



imaginary part of the action, expression (|13[) becomes 

f S} dgi 
S= / Pi(s)-^-(s)ds-£:ReT 

+ i^J% 1 ( S )^( S )ds + E(-ImT^j . (51) 

Then 



S e (E) = S (E)-ET 



(52) 
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with T = T and 



S {E) =w,S z (E)+iS c {E)/2 

= 2w t [ y / 2(E + -/cos <?) d <? 
Jo 



2i 



rgch(£;/ 7 ) 



y / 2(E-jchq) dq . (53) 



All these expressions can be written in terms of the com- 
plete elliptic integrals [37j (defined for |it| < 1) 



JC(u) ^ / -= 
Jo VI 



d.r 



9 • 2 

u z sin a; 



/.?r/2 

= / \fl — u 2 sin 2 x dx . 
Jo 



Namely, 



S Z (E) = A^2(E + 1 )£ 



S C (E) = 8^2(^ + 7) 

Then we obtain 
A, 

with 




(54a) 
(54b) 

(55a) 



-E 



2h 



3 -S c (£„)/(2ft) 



T t (E) = 



T x (E n ) 
2V2 



' E — j 
£ + 7 

(55b) 



(56) 



dS r 



dE 



K 



27 
£ + 7 



(57) 



The energies of the highly excited states are approxi- 
mately given by the free rotations: E n ~ n 2 h 2 /2 S> 7. 
Then T z {E n ) ~ 2n/(nh). The asymptotic expansion 
of <S C (.E) for large leads to 



A, 



ft 2 V2S 



it h 
1 /e 



n e 2(ln(7/£„)+2-31n2)\/2~E^//i 



7rn' 



4n-l 



(58a) 
(58b) 



The last expression corresponds exactly to equa- 
tion (3.44) of |35l ] obtained with standard uniform semi- 
classical analysis. We see on figure [5] that (1B"6|) is a very 
good approximation even when the energies E n get close 
to 7 the energy of the separatrix. 



VI. RESONANT TUNNELLING AND 
FABRY-PEROT EFFECT 

We are now ready the see how formula (|4T)|) allows us 
to reproduce the resonant tunnelling between two wells 
related by parity when the potential in (TTJ) has also a 



n = 1 



n = 9 




12345678 



b) 




6 7 8 



Figure 5: For the pendulum with 7 = 1, we have plotted 
in a) the exact value of ln(Ai? n ) (see eq. (|47[0 versus l/h 
for n = 1, ...9 (black lines). The dots are given by the 
semiclassical approximation ln(A n ) (see eq. (|58bp ) up to the 
maximum value of l/h (vertical dotted lines) where be- 
come lower that the separatrix energy 7. The semiclassical 
limit is obtained when the states become more and more lo- 
calised in the region of phase space corresponding to rota- 
tions: for a fixed n this requires to increase E n ~ n 2 h 2 /2 that 
is to increase h; or for a fixed energy E > 7, we must in- 
crease n oc V2E/h. To control the validity of the prefactors, 
we can check in b) that ln(A_E rl /A n ) approaches zero in the 
semiclassical limit. 



deeper central well (Fig. [5] a) . The minimum of the right 
and left wells is fixed at zero, the minimum of the central 
well is V(0) = Knin < and the local maximum between 
the wells is denoted V max . When < E < V max , we 
will denote by q r (E) < q' r (E) < q"(E) the three posi- 
tive solutions of V(q) = E. As explained in section (IIV[) . 
we will try to construct appropriate time paths [t] to 
exhibit complex trajectories with purely real q that 
connect the two symmetric tori from t to I at some 
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Figure 6: (colour online) To describe resonant tunnelling be- 
tween two symmetric wells centered at ±a and separated by 
a third central well (a), we use orbits made of five primitive 
trajectories with a time path given in (d). We still maintain 
the position q to be real and can use a three dimensional sec- 
tion of phase-space (see text). Here we have chosen V given 
by ([65]) with a = 7/4 and 6=1/2. 



energy < E < V max . These tori are delimited by the 
two turning points ±p' — (0,±q' r ) and ±p" = (0, ±g"). 
What is new of course is the existence of a central real 
torus m delimited by p — (0, q r ) and —p. Using the 
three dimensional, representation of the section Im q = 
of phase space (Fig. [5] c and e), we see that the or- 
bits o must be a series of concatenation of five trajec- 
tories connecting at the turning points ±p, ±p' . First 
we start with one portion living on r with dt/ds > and 
Imp = 0, then connect at p' to a trajectory with Rep = 
where idi/ds > 0. It follows the energy curve c whose 
equation is (Imp) 2 /2 — V(q) = —E. Then o can con- 
nect at p to a real trajectory on m with Imp = 
and dt/ds > 0, then can cross — c from — p to —p' be- 
fore reaching [ = — r. The corresponding time path will 
necessary have at least two steps (Fig.[6]d) each of them 
having a height which is a half- integer multiple of T c , the 
real period of the primitive periodic orbit c. Among all 
the possible o's, we will keep only the exponentially dom- 



inant contributions, when o remains as shortly as possible 
with complex p. Then, for such orbits to exist, we must 
choose T of the form 



T = t + w x T z (E) 



\)T m[ E)-xT A E) (.-,!» 



involving the periods of the primitive orbits and the 
corresponding winding numbers w Xl w m which are non- 
negative integers. r denotes a positive fraction of 
time smaller than T Z (E). The base points for the two 
branches (3 = ± defined by Rep ^ on t coincide with 
the turning point b± = q' r , the action S (bpr,bp,E n ) and 
the index p are independent of the choice of the branch 
where o starts: 

S (bp,,b p ,E n ) =w t S v (E)+(w m + ^j S m (E) + iS c (E) . 

(60) 

(61) 

(62a) 
(62b) 
(62c) 



p = w v + w m + 3 



Explicitly, we have 



S t (E) c = f 2 



S c (E) d ^2 



1'AE) 
q r {E) 



q r (E) 



2{E-V{q)) dq- 
2(E - V(q)) dq ■ 



2(V(q) -E)dq; 



and the corresponding periods T r , T m , T c are obtained 
by deriving with respect to E. 

As we have seen, formula (|40[) will provide an approx- 
imation of the exact splitting AE n that becomes better, 
the better the condition ReT(E n ) 3> ImT(£ n ) is sat- 
isfied. Not only this condition render the precise value 
of t irrelevant, it also requires large w m and/or u; r , espe- 
cially if we work with E n > for which T c (E n ) > T m {E) 
and T c (E n ) T t (E). For a given pair of w v and u> m , 
there are 4(u> r + 1) topological classes of orbits corre- 
sponding to two possible initial branches, two possible 
final branches and r = 0, . . . , w t possible windings on r 
for w t — r windings on I. Then all different w t and w m , 
such that (j5U)) holds, give a contribution to (T4"0"|) : If we 
define R(T) = Re T — r — T m /2 ~ Re T then, 



A»(T) 



h^o T 



^ e -s c (E n )/h 



X 



2^ K+i)c 

{wz,w m } pos. int. such that 
w t T t (B„)+u; m T TO (£;„)=iJ(T) 



iw v [S v (E n )/h-w]+iw m [S m (E n )/H-Tr] 



(63) 



We immediately see the resonance at work since the sum 
reaches a maximum when both v x =S v (E n ) / \2ttH) - 1/2 
and u m =S m (E n )/(2nh) — 1/2 are integers: the energy of 
a state mainly localised in the central well becomes nearly 
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obtain 




Figure 7: (colour online) For V given by (|65|l with a = 7/4 
and b = 1/2, we plot ln|A_E n | (thick black line) for the low- 
est doublet in the symmetric wells ~ hw/2 with lo = 
2a\j2{a? — b 2 ) ~ 6.30. The (blue) dots correspond to estima- 
tion (jXT3|) with T = \T\e- ie = (K + \/2)T m {E n ) -iT c (E n ) in 
the limit K — > oo (9 — >■ 0). The upper thin (red) curve pro- 
vides In | AE n / A n \. The inset provides a magnification of the 
third spike around 1/h = 6.6 for which T m (E n ) = 2T t (E n ) ~ 
1.58 and T c (E n ) ~ 1.50. Here, we have also plotted In |A„(T)| 
for if = 100 (thin orange (light gray), 9 ~ 0.01) K = 30 (thin 
green ~ 0.03), K = 5 (dashed thin magenta 9 ~ 0.17). 



degenerate (up to ft 2 terms) with the energy doublet in 
the lateral wells. Then the contributions of the repeti- 
tions of m interfere constructi vely like the optical rays in 
a Fabry-Perot interferometer [2ll sees. 12.14-12.17]. To 
estimate the sum in the right hand side of (|63[) , let us take 
a rational approximation of the ratio T m /T r , namely 



m 
r 



(64) 



with to and r being coprimes positive integers. For the 
polynomial potential 



(65) 



with a > b > 0, the argument presented in [38[ can be 
generalised to show that ([64]) is actually exact with r = 1 
and to = 2 for any energy < E < V max . If K denotes 
the integer part of R(T)/(rT m ), we can compute and ap- 
proximate for K ^> 1 the right hand side of (1551) and 



|A„(T)| 



-S c {E n )/h 



h->0 T v | sin (7r(m!/ t — rv m )^ | 



(66) 



Fig.[7]shows that this latter expression provides a good 
approximation for AE n even in the immediate neigh- 
bourhood of a resonance where estimation (fT0|) is not 
justified any more. If we had continued working with a 
finite K, the sum (|63[) would have involved a finite num- 
ber of terms and the singularities due to the vanishing 
denominators in (1661) would have been smoothed down 
(inset in Fig. [7]). In other words, for a fixed \T\, rotating 
down T in the lower half plane ImT < 0, destroys very 
quickly (0 ~ 0.2) the large resonant fluctuations of tun- 
nelling. This effect has already been shown in the case 
of a kicked system [28]. 



VII. ESCAPE RATES 

So far we have focused our analysis of tunnelling in 
bounded systems only, but the philosophy we presented 
here can be extended to more general situations. For ins- 
tance, let us show how we can compute the escape rate 
from a metastable state localised in a confining poten- 
tial V whose shape has the form given in Fig. [S]a). The 
potential has a local minimum at q — (say Vmm = 0) 
and an energy barrier for < q < g ma x whose height 
is V max . For q > q max , the potential remains non-positive 
and therefore in real phase-space (Rep, Reg), V defines 
around the origin an island of stability made of tori with 
positive energy. One state whose Husimi distribution is 
initially localised in the island, say a quasi-mode |$ n ) at 
energy < E n < V max , will progressively decay outside 
the well. The decay rate T n is then defined from the 
overlap: 



_ -r n T/2-iE„T/h 



(67) 



(* n \U(T)\$ n )=e 
If we choose a complex T such that 

|r„T| « i , 

then we obtain a trace formula for T n with the help of 
the projector-like operator Tl n 



(68) 



r„ ~ 



ImT 



■ Im 



(e^ T 'Hr(n n U{T))) (69) 



which allows an explicit semiclassical expansion in terms 
of classical solutions with a complex time path. The dom- 
inant contributions will be provided by periodic orbits in 
time T with real q starting on the torus at energy E n , 
then, while idt/ds > 0, going forth outside the well be- 
fore coming back to its initial starting point (Fig. [8] c, 
d)). Then, up to a dimensionless factor f n of order one, 
we get 



fn 



T c (E n ) 



-S c (B„)/ft 



(70) 
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Figure 8: (colour online) The escape rate at energy E from 
an island created by a confining potential V whose shape is 
shown in a) can be semiclassically computed in terms of com- 
plex classical orbits where one canonical coordinate, say q, is 
kept real. In phase-space, this trajectory appears to be a con- 
catenation of two types of curves that join on the (Reg)-axis: 

1) a trajectory that lies in the phase-space plane (Reg, Rep) 
with real variations of t at energy E with the potential V and 

2) a trajectory that lies in the phase-space plane (Reg, Imp) 
with imaginary variations of t at energy — E with the poten- 
tial — V shown in b). In c) a family of constant energy curves 
is shown (horizontal green for the first type, vertical red for 
the second type). In d) for a given energy, we show how three 
curves (, c, r glue together at the turning points. 



with 



S C (E) d = 2 



q T (E) 



2(V{q) - E) dq (71) 



where < q r {E) < q' r {E) are the two right turning points 
at energy E. Here we have supposed that the action be- 
tween the two left turning points is larger than ([71]) and 
gives a subdominant contribution. If the potential is sym- 
metric, two symmetric complex orbits would contribute 



with the same weight and therefore the escape rate should 
twice as large. In the case of an island with sharp bound- 
aries around g max , S C (E) which is the area enclosed by 
the primitive orbit c (Fig. [5] d) is mainly given by the 
portion of c where we can keep an harmonic approxima- 
tion for the potential: V(q) ~ w 2 q 2 /2, q r (E) ~ y/2E/uj, 
q' r (E) ~ \J A/ (ttu}) where A is the area of the island in 
the real phase space. Then, with a =ujA/(2ttE), we have 



2E 



S C (E)~ — [vWa^-lnfVa + Va 3 !)] • (72) 

id 



Inserting this expression in with E n ~ (n + l/2)fiw, 
we exactly recover the expression (5) used in [39j j with an 
elegant and simple interpretation. The chaotic sea that 
surrounds the integrable island in the mixed system con- 
sidered by Backer et al. acts as a sharp effective potential 
barrier as the one draw in picture [8] our complex trajec- 
tory that allows to escape from the regular region has its 
main features governed by the integrable (and even har- 
monic) approximation of the dynamics about the island, 
following precisely the general philosophy of [3!| and [4(| ■ 
This computation of the "direct" tunnelling (by opposi- 
tion to resonant tunnelling where the model of a pure 
quadratic kinetic energy fails) can also be reproduced 
within the standard one-dimensional JWKB theory used 
for computing transmission coefficients. 

Here again, we can check easily that the tradition- 
nal instanton method is included in our approach : The 
regime dlT — > selects the instanton solution and pro- 
vides the escape rate from the equilibrium point as given 
by equation (2.47) of 2, chap. 7]. 



VIII. CONCLUSIONS 

The explicit semiclassical expansions of trace formu- 
lae for tunnelling splittings (or escape rates) in terms 
of classical orbits constructed with complex-time paths 
provide an interesting alternative approach to Herring 
formulae essentially because they do not require to an- 
alytically continue the wave functions in the complex 
plane. In multidimensional tunnelling and/or for a non- 
autonomous Hamiltonian system, the generic lack of con- 
stants of motion isolates the stable islands (if any) from 
each other by chaotic seas. The analytic continuation 
of the kam tori that build the islands is prevented by 
the existence of natural boundaries (see for instance the 
recent discussion in [30 . part II, appendix B] and ref- 
erences therein). The approach we have presented here 
seems to circumvent these difficulties but, of course, the 
problem of how to select, with an appropriate [t], the rel- 
evant trajectories among an a priori exponentially grow- 
ing number of classical complex solutions remains open. 
We expect the tunnelling splittings between two states 
at energy around E localised in two symmetric islands to 
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be approximated by the expansion of the form 



I T 



[Pl,P2,---] 



J 1 w i lSi(E n )/h-ir] - J2f = i w 'j[S'j(En)/h] I 



(73) 



where the sum runs over all possible sequences of turning 
points [pi , p2 , ■ ■ ■ ] at energy E such that we can choose 
a complex-time path that leads to a trajectory, made of 
primitive orbits, that connects in time T the two (real) 
tori. The winding {wi} and the actions {Si} (resp. {w'j} 
and {Sj}) refer to the primitive orbits obtained when the 
variations dt/ds are purely real (resp. purely imaginary). 
For dimensions larger that one, the dimensionless pref ac- 
tor / may appear as a power law in h. Here, inspired by 
the study in section PvTl we can qualitatively see how the 
constructive interferences between repeated paths emerge 
in a speckle- like forest because of the presence of resonan- 
ces. As shown in [28| . a progressive complex rotation of 
time provides a natural way to select the main resonance 
effects. If we want to expand the splittings (or the es- 
cape rates) according to elementary process as proposed 
in [40| , our approach offers a promising tool to interpret 
and compute semiclassically all of the ingredients of such 
an expansion. 
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Appendix A 

In this appendix, we derive the dominant contribu- 
tions (flU)) to the traces involved in ©. Both cases 
tr([/(T)) and tr(S U(Ty) will be treated simultaneously 
by defining the sign -q to be +1 in the first case and —1 
in the second case. The semiclassical arguments under- 
pinning the derivation are relatively standard and may 
be found in one way or another in the literature. For 
instance, the contribution (|A13|) for r\ = 1 can be found 
in I4ll. e g. (2.12)] within a more restricted context (see 
also [42, chap. 8]). Nevertheless we found it useful to pro- 
vide all the steps in the precise context of this work, not 
only to render the presentation self-contained, but also 
because we are working in the time domain with general 
Hamiltonians that have not necessarily the form (fT]). 



Given a [t], for any classical phase-space path o, we can 
consider the final coordinates (pf,qf) at time t(sf) = T 
as smooth functions of the initial coordinates (pi,Qi) at 
time t(si) = 0. The monodromy matrix M is defined 
as the differential of these functions: to first order, an 
initial small perturbation implies the final perturbation 



f5p f 
\Sq f 



M . n M 0i ia 

M oi M„ 22 



(Al) 



The sub-matrices can be expressed from the sec- 
ond derivatives of S (qf ,qi\T) by differentiating equa- 
tions CEBD: 



i^pqjqi^oj 



M„ A1 = -8$ fqf S l 



Mo, 12 — dq fq .S a — dg fqf S (dq fq .S c 



M ,21 
M , 22 



ci 2 S" 

If H ° 

ci 2 <5 

If 1i 



8 2 • 

' 



provided that d 2 fq .S is invertible. Equations 
can be inverted in 



d. 



if if 

while (|A2b ) provides 



q [qi So — -(M 0t 2\) ; 

d im S o = (M 0)2 i)- 1 M ,22 
dl qf S =M , 11 (M 0i21 )- 1 



(A2a) 
dl qi S ; (A2b) 
(A2c) 
(A2d) 
2k,c,d) 

(A3a) 
(A3b) 
(A3c) 



M , 12 = -(M^i)- 1 + M , 11 (M 0)21 )- 1 M 0! 22 (A4) 

which is nothing but the expression that det M — 1 once 
we use the identity 



det 



A B 
C D 



detiCAC^D - CB) (A5) 



where (A, B, C, D) are square matrices of the same size 
and C is invertible. 

Coming back to the oscillating integral p7|) . if the in- 
tegration path can be deformed in order to pass through 
an isolated critical point q c of q i— > S (r)q, q,T), we will 
have a contribution of the form 



'det 



( d 2 S 



\dq l dq f 



\ 



det - 



d 2 (S ( m ,q;T)y 
dqdq 



(A6) 



where all the functions are evaluated at {rjq c ,q c ,T). 
Physically q c is interpreted as the initial position of the 
phase space paths o such that (pf, qf) = (ijpi, r/qi). With 
the help of (|A3[) , the prefactor of the exponential can be 

written as [det (M ,22 + M ,2iM ,ii(M ,2i)" 1 -2t])]~ 1/2 
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up to a global sign. Using 
tion (|A6[) can be written as 

(-1)". 



5J) again, the contribu- 



V(-?7) D det(l-77Af ) 



iS„/h 



(A7) 



where (— 1)"° fixes the sign of the square root and D de- 
motes the number of degrees of freedom. Recall that in 
the definition of the path- integrals (fT2|) . a time-slice r is 
always implicit. Let us keep for a moment an explicit 
discretization, for instance with T being an integer mul- 
tiple of r, o referring to a discrete set of points, Hamil- 
ton's equations (|T1| being discretized into a phase-space 
map, path integral (TT3")) being turned into a discrete (Rie- 
mann) sum, etc. Then, the contribution of each of the 
points of a given o, such that (f>/,5/) = (vPiiVQi)] is the 
same and remains given by (|A7[) ; but this is correct only 
if, for a given r, h is small enough, because when per- 
forming the steepest-descent method on (fT?]) , one must 
be able to split the oscillating integral into separate con- 
tributions coming from two distinct points of o. In the 
continuous time limit, i.e. when the limit r — > is taken 
before the semiclassical limit H — > (see [??]), the or- 
bits o of non-zero length appear as a one-dimensional 
continuum none of whose points can be considered sep- 
arately anymore. The only isolated critical points q c are 
given by the position of some equilibrium points. Under 
the symmetric condition (qi,Pi) = (~Pi, —Qi), only the 
origin must be examined. Linearising the Hamiltonian 
flow about a non-degenerated fixed point e = (p e ,q e ) 
leads to a monodromy matrix whose eigenvalues can be 
collected by pairs exp(±A CiQ T) where {X^ a }a=i,...o are 
the Lyapunov exponents. Then (|A7|) becomes 



(A8) 



(-1)"' e -iH{Pe,qe)T/K 

nLi(e A - T/2 -^-^ T / 2 ) 



with a possible adjustment of the sign. For a generic 
choice of T, the denominator does not vanish. 

For a non-zero length path o, the critical g's are de- 
generate along the trajectory; for a system with several 
degrees of freedom, one must treat separately the (gaus- 
sian) integrals on the transverse coordinates q± along 
which S (rjq,q,T) varies (quadratically) from the longi- 
tudinal coordinates q\\ along which S (r)q, q, T) is cons- 
tant. Of course, the dimensions of q± and q\\ depend 
crucially on the presence of KAM tori. However, multi- 
dimensional tunnelling is beyond the scope of this paper 
and, the quantitative studies presented here concern one- 
dimensional systems only. The contribution to the trace 
of such a path is then, up to a global sign, 



dq* 



d 2 S 



dqidq f 



(A9) 



The conservation of energy along o, H(p, q) = E, implic- 
itly defines a function p(q, E) in the neighbourhood of any 
point where d p H ^ 0. Globally along the trajectory o, 



we may encounter several possible branches pp(q, E) for 
the graph of these functions (the two possible signs of 
a square root when the Hamiltonian has the form ([T])) 
which become singular but pairwise connect smoothly 
at the turning points, defined by d p H = 0. A relation 
between E, T, qf and qi can be obtained by integrat- 
ing (dt/ds) ds and using (I14b[) : 



T = T o(qf,Qi;E) 



dcf 



dq/ds 



d p H(p(q(s),E),q(s)j 

(A10) 

Everywhere but at the turning points, the value of s 
dictates the choice of the branch used for the inte- 
grand. This relation implicitly defines E (qf,qi;T). 
The usual expression for the derivative of implicit 
functions leads to the relations: dEP — l/d p H and 
d qi E = -d qi T /d E T = d T E /d p H(p i ,q i ). If we dif- 
ferentiate (|18bl) with respect to when pf is given 
by p(qf,E (qf,qi,T)), we obtain 



d 2 S 



dqidq } 



1 



1 



dE 



d p H(p f ,q f ) d p H(pi,qi) dT 



{if,n;T) 
(All) 

With property ((3]), we get d p H(r]p,T)q) — r]d p H(p,q) 
and then, by differentiating (|A10p . we have 
d q (T (riq,q, E)) = 0. Therefore, the energy E of 
the path o depends only on T, not on its starting 
point q. The square root of dTE a (rjq,q,T) = dE /dT 
can be got out from the integral in (|A9[) . When adding 
the contribution of each path whose starting point lie on 
the branch /3, we obtain 



Tb{E) 



clef 



dq 



d p H[pfi(q 7 E),q 



(A12) 



This is not exactly the right hand side of (|A10[) because 
the domain of integration in (IA12[) is the domain of the 
branch where the starting point of o lives. Each branch /3 
is delimited by two turning points and Tp is the time 
spent to go from one point to the other. 

The integral (IA9|) involves all the possible starting 
points for a trajectory o and therefore we must add all 
the branches that patchwork smoothly in phase-space to 
form the geometrical set of points crossed by o. Referring 
to the purely geometrical quantities (i.e. independent of 
the choice of the parametrization), we have the contribu- 
tion 



i„ V 



v 7 — 2r]mh V dT 



dE, 



o e iS /n 



(A13) 



only if T = T: for rj = +1 the path o is a periodic orbit 
and, for rj = — 1, the path o is half a symmetric perio- 
dic orbit (the whole periodic orbit being of period equal 
to 2T) . The sum concerns all the geometrical branches /3 
crossed by o (even if o passes several times by the same 
points, each branch is only counted once). As before, the 
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conversion of a product of two square roots of complex 
numbers to the square root of the product may introduce 
a sign that can be absorbed in the definition of the 
exact computation of the index fj, is difficult but since 
it may change at the bifurcation points only, where the 
scmiclassical approximation fails, it is sufficient to know 
that it depends on the nature and the number of the 
turning points encountered on o. Therefore it is an addi- 
tive quantity when several primitive orbits are repeated 
or concatenated together. 



Appendix B 

In this appendix we explain how to obtain the asymp- 
totic expansions (|33| as E — » + . 

First consider S x and split it in two parts S+ + S- 
whcrc 



S + d ^ c 



2J2(E-V(q))dq; (Bla) 



clef 



SL = / 2J2(E-V(q))dq . (Bib) 

Jq r {E) V 

Setting e = a — q r {E), rewrite S- as 

S_ = 2v / 2e / ^V(a-e) - V(a - sejds, (B2) 
Jo 

Now expand the integrand as a power series in e up to the 
fourth order, compute the integrals that appear in each 
coefficient and insert the expansion of e in E obtained 
from the implicit equation V(a — e) = E: 



V2 m V^(a) „ 
e = — VE+ E 



4^ (4) (»)-5^ (3) («)) : 



36w 7 



E 3/2 + 0(E 2 



(B3) 



Proceed in an analogous way for the computation of the 
first three terms of the asymptotic expansion in \[E for 
S + . When summing S+ and S-, expression (|33b[) is ob- 
tained with (1331). 



The expansion of S c is more subtle since it is not dif- 
ferentiable at E = 0. Its derivative is given by 



dS c 
~dE 

def 



J L(s,e)ds (B4) 



where we denote e = a — q r (E) and define 
def 2y/2 



L{s,e) 



y /V(a-s)-V(a-e) 



(B5) 



The function L{s, e) is not continuous but we can extract 
the discontinuous part from 



L(s,e) 




(-l)"2t/(™)(a) s 7 



by expanding the last factor: 

4 2V^(a) 



L(s,e) 



(B6) 

+ M(s, e) 
(B7) 

where now M(s, e) is a continuous function of its two 
variables. Then, a standard theorem in analysis as- 
sures that e i — y J a M(s, e) ds is continuous and its limit 
when e — > is 



LOWS 2 



3w 3 ( fl 2 _ e 2)3/i 



M(s,0) ds 



L(s,0) 

o V ws 

4A 2aV { - 3 \a) 

id 3U! 3 



3w 3 



ds 



(B8) 



with A given by ([34]) . The two other integrals obtained 
by inserting (|B7[) in the right hand side of (|B4[) can be 
computed exactly and expanded as e — > + up to or- 
der o(l). Then, inserting (IB3|) . we obtain 



(B9) 




Its integration leads directly to (|33ap . 



Appendix C 

The quasi- mode |$ n ) ==(|<?^) + |0^))/\/2 is localised 
on one torus at energy £7+ ~ E~. Standard JWKB tech- 
niques [H, 13] provide a semiclassical approximation to 
its wave function: 



.4 



n,0 



J2fj T p( E n) Jd p H(pp(q,E n ),q) 



x exp \ij pp(x,E n )dx/hj (CI) 

(/3 labels the possible several branches of the torus, A n ,/3 
are dimensionless coefficients of unit modulus, bp is a 
base point of the branch j3 and Tp(E n ) the character- 
istic time (|A12j) spent on the branch j3). Within the 
scmiclassical approximation, it is, therefore, consistent to 
construct II n (q' , q) by substituting (|C1|) in the matrix el- 
ements of the projector operator |</>+) (<^+| + \cj>~) (<f>~\ = 
|$„) ($ n | + S |$„) ($ n | S. From the integral O 

2 J dqdq' $r (<7) ($r (9'))*G( W ' !g ;T) , (C2) 
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when we insert the semiclassical expressions (|16l) , we ob- 
tain a sum of integrals of the form 



dq dq 



p,3{x,E n )dx / p '(x,E n )dx 

b p e n Jb fl , 



^d p H(pp(q,E n ),q) yj ' d p H{pp,{r)q' ,E n ),q') 

e iS (nq',q;T)/H _ /Qg\ 



y %6% 

The stationary conditions 

Pp{q,E n ) = - 
r)Pf3>(q',E n ) = 



(m',q\T) 



d qi S (r)q',q;T); (C4a) 
d qf S (m',q;T) (C4b) 



select the classical trajectories o with energy E n that go 
from (pi,qi) = (pp(q,E n ),q) at t(si) = to (j)f,qf) = 
(rippi(q' ,E n ),r)q') at time i(s/) = T. Then the value of 
the exponent 

S 0lP ,p>(E n ,T) = -E n T + S (bp,,bp,E n ) (C5) 

depends only on the branches where the starting and 
ending points lie and not on the precise location of these 
points on the branches. Since q and q' correspond to 
the same torus, such a trajectory must connect the two 
symmetric tori for 77 = — 1. At a given E n and T, for 
a fixed q = qi on the branch /3, qf and /3' are uniquely 
given and we can make the stationary phase approxi- 
mation for the integral on q' . Then if we insert (|A10[) 
at energy E n into (|C4I) and differentiate it with respect 
to q or q' , we obtain some identities that, with (jAllI) . 
allow us to simplify the combination of the prefactors 
and the remaining integral in q turns out to be precisely 
of the form of the right hand side of (IA12[) . A priori, 
the domain of integration is included in the domain of 
the branch (3 but is not necessarily equal to it because 
when sliding the starting point on the whole branch f}, 
the endpoint may cross a turning point and correspond 
to a jump of /?'. However, we obtain characteristic times 



that depend only on the geometry of the orbit, not on the 
number of times the considered branch may be repeated 
as s goes from s» to s/. As discussed in the case of the 
double well, if there exist different topological classes of 0, 
each of them being characterised by an ordered sequence 
of turning points [pi , p2, ■ • ■ ] , we must add such contri- 
butions. Then, using directly tt(n n U{T)) ~ 2e-' lE ™ T / h , 
we have proven that 

A «( T ) ~ % E (-ir^(£„,T)e i5 °(V^^)A 
ft->0 1 ' — ' 

[pupa,...] 

where the sum runs over all the sequences of turning 
points on the section s at energy E n where one canonical 
variable is maintained real. There must exist for such 
a sequence, one half symmetric orbit starting on the 
branch /3 of the torus at energy E n , crossing successively 
all the sequences [pi, p%, . . . ] and ending on the branch /3' 
at time T. The dimensionless coefficients Ap^>(E n ,T) 
have a ft-independent modulus of order one and depend 
only on the geometrical properties of the branches. If 
some parts of the trajectory are repeated, their repeti- 
tion numbers do not appear in A^^i(E n ,T) but only in 
the cumulative quantities: the index /i and the action 
S (bff',bff,E n ) given by 

f s f dq 
S (bp>,bp,E n ) = J p{s)—ds 

pp(x,E n )dx- p0>(x,E n )dx. (C7) 

In the case of two branches, the computation of the coef- 
ficient can be done exactly using the appropriate choice 
of phase conventions for the base points bp and A n .p; we 
obtain: 



A 0t0 ,(E n ,T) = -1/2. 



(C8) 



We illustrate in the main body of this article, how to 
compute the sum in the right hand side of (|CJ6|) . 
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